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We study spin imbalance effects on the Larkin-Ovchinikov-Fulde-Ferrel (LOFF) state relevant for 
superconductors under a strong magnetic field and spin polarized ultracold Fermi gas. We obtain 
the exact solution for the condensates with arbitrary spin imbalance and the fermion spectrum 
perturbatively in the presence of small spin imbalance. We also obtain fermion zero mode exactly 
without perturbation theory. 
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I. INTRODUCTION 

The exotic superconducting state called Larkin- 
Ovchinikov-Fulde-Ferrel (LOFF) state has been proposed 
to arise in superconductors under a strong magnetic 
fieldjii^ The LOFF state has a spatially varying order pa- 
rameter associated with Cooper pairs with finite center- 
of-mass momentum. If a strong magnetic field induces 
spin polarization, such Cooper pairs are considered to 
form between electrons with different Fermi momenta. 
The LOFF state is also relevant for the physics of color 
superconductivity where quarks with different masses 
form pairsi^ This state was not observed for over 40 years 
since its proposal, in spite of tremendous efforts. In the 
last couple of years, there have been several claims of 
its observation in heavy fermion materials^ and organic 
superconductors^, but direct confirmation is yet to be 
given (see Ref. IH for a review) . 

Recent developments of research in cold atomic Fermi 
gases have renewed interest in the LOFF state (see Ref. 
for a review) . In two component Fermi gases consisting of 
atoms in two different hyperfine states, (pseudo-)spin po- 
larization can be controlled by changing the populations 
of the two components »i Furthermore, atomic interaction 
can be tuned in this system by using the Feshbach reso- 
nance which allows one to explore the interesting BCS- 
BEC crossover physics. Thus, a spin polarized Fermi gas 
is an ideal system for realizing and exploring the LOFF 
state. In Fermi gases in a toroidal trap, it has been shown 
that a new state called angular LOFF state is possible in 
which the rotational symmetry is spontaneously broken;^ 
instead of the translational symmetry for the usual LOFF 
state. Recently, observation of spin polarized superfluid 
state was reported^ and it is expected that the LOFF 
state has been achieved in this experiment. However, di- 
rect observation of its oscillating order parameter is still 
lacking. 

The Bogoliubov-de Gennes (BdG) equation has been 



widely employed to study the LOFF state. Machida and 
Nakanishil- derived the self-consistent LOFF state so- 
lution of the ID BdG equation making use of the an- 
alytical solutions of the ID Peierls problemi^ir— How- 
ever, they assume that electrons with up and down spins 
have the same Fermi velocities {vf^ = v-p^), so that their 
LOFF solutions are valid only when spin polarization 
is small. This assumption is appropriate to supercon- 
ducting states, because in ordinary superconductors the 
splitting of the Fermi surfaces is of an order of the pair 
potential at the Pauli limit, which is much smaller than 
the radius of the Fermi surfaces.^" On the other hand, 
the Fermi surface mismatch is in general not small for 
cold atomic Fermi gases^ and therefore we have to take 
into account large spin polarization. The spin imbalance 
effect was previously studied in the Peierls problemJ^ 

Recently a new approach for solving the BdG equation 
has been proposed by Ba§ar and Dunne. They derived 
the non-linear Schrodinger equation (NLSE) for the order 
parameter A with a suitable ansatz for Gor'kov Green's 
function. Since the derived NLSE is a closed equation 
for the order parameter A(a;), this enables one to avoid 
the self-consistent calculation of the coupled equations 
of the BdG equation and the gap equation. Using this 
approach, they found a new self-consistent solution for 
a complex kink crystal, which includes all previously 
known solutions as special cases, such as the solutions 
of the LOFF state (real kink crystal)^ and Shei's com- 
plex (twisted) kink^^. This new approach and the com- 
plex kink crystal solution have been further developed— 
for the massless Gross-Neveu modeli^ and the Nambu- 
Jona-Lasinio model in 1 -|- 1 dimensionJS However, this 
approach has not been extended to spin polarized system. 

In this paper, we investigate the self-consistent solu- 
tions of the BdG equation for spin imbalanced Fermi con- 
densates. We extend the approach developed by Ba§ar 
and Dunne to obtain the analytic solutions for the LOFF 
state in which the order parameter exhibits spatial os- 
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cillations. In contrast to the solutions of Machida and 
Nakanishi, we take into account the difference in the 
Fermi velocities (uf^ ^ wp^) and derive the exact solu- 
tions for the condensate A(a::) which are valid for any spin 
polarizations. We show how the effect of large spin po- 
larization changes the form of the non-linear Schrodinger 
equation for the order parameter. We also develop the 
perturbation theory for the BdG equation in the pres- 
ence of small spin polarization, and obtain the fermion 
zero mode which is exact for arbitrary spin polarization. 
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where v-pa = kp^/m is the Fermi velocity and A ~ 
e-*('=ft+fcF4-)^Ao. When vp^ = vpi, the LOFF state solu- 
tion of Eq. dH) has been derived in Refs. HoIITtI . 



B. The Ba§ar-Dunne formalism 



II. NONLINEAR SCHRODINGER EQUATION 
FOR ORDER PARAMETER 



In this section, we derive the nonlinear Schrodinger 
equation for the order parameter A in the presence of 
spin polarization. 



A. The Bogolibov-de Gennes equation 

We consider a gas of fermions with spin up and down in 
quasi-one dimension under a magnetic field. If fermions 
with different spins interact attractively, the system un- 
dergoes a superconducting (superfluid) transition at low 
temperature. Although the mean field approximation is 
not valid in strict one dimension, since we assume a quasi- 
one dimensional system relevant for experiments, the sys- 
tem can be described by the mean-field BdG equation^^ 
(we set h= 1) 
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where (T(=t,|) stands for the spin and M is the mass of 
the fermion. The energy difference due to the Zeeman 
splitting is included in the chemical potential for each 
spin state fia (cr =t,i). This model is indeed applica- 
ble to an imbalanced cold Fermi gas;^ Throughout this 
paper, we restrict ourselves at T = 0. In this case, the 
order parameter Ao(a;) satisfies the gap equation 



Ao(a;) = -2g^ ^ u„(a;)v„(a;) 



(3) 



where g is the attractive interaction between fermions 
with different spins and n is the index for eigenstates. 

If the attractive interaction is small compared with the 
Fermi energy epo- — ^J■cr , fermions near the Fermi surfaces 
form Cooper pairs. If we assume uo{x) — e^''^^^u{x) 
and vo{x) ~ e~'^'^'^^'^v{x) (fepo- is the Fermi momentum 
/cpo- = \/2Msp^), u{x) and v{x) vary much slower than 
1 /fcpo-. Neglecting the second derivative term of u{x) and 
v{x) (the Andreev approximation^), the BdG equation 



In Ref. [13, the so-called nonlinear Schrodinger equa- 
tion (NLSE) for A(a::) has been derived for the case of 
v-p-f — vpi, through the analysis of Gor'kov Green's func- 
tion. This is a convenient way to solve the BdG equa- 
tion. Since the derived NLSE is a closed equation for 
A(a;), this enables one to avoid the self-consistent calcu- 
lation of the coupled equations of the BdG equation and 
the gap equation. We extend this analysis to the case of 

First, we derive the Gor'kov Green's function that sat- 
isfies 



{H-E)G{x,y;E)^S{x^y), 



where 
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-ivF^^ A{x) 
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The Gor'kov Green's function can be constructed from 
two independent solutions ip{x) and (j){x) of Eq. ^ as2^ 
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V- 
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where W = iip'^a2(f> is a Wronskian. It is easy to show 
that the Eq. ^ satisfies ([5]). 

On the other hand, the diagonal resolvent is defined 

by 



R{x;E) 



1 



H -E 



(9) 



Indeed, Eq. (|9]) includes all spectral information for 
fermions in the presence of A(a;), such as the single- 
particle spectral function 

p{E) = ilm J dxTrR{x;E + i5). (10) 

From Eq. (O, R{x; E) can be obtained as the coinci- 
dent limit of Gor'kov Green's function: 

R{x;E)^ lim l[G{x,x + S;E) + G{x + S,x;E)]. (11) 
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We can show that R{x) satisfies the following conditions: 



TrR{x; E)a3 = 
det R{x;E) = 



where 
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In addition to the above conditions, from the definition 
Eq. (|9]), the resolvent must satisfy the Hermiticity con- 
dition: 



R\x]E) = R{x;E). 
By a straight forward calculation, we obtain 



(15) 



d^Rix; E)a3 
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,R{x;E)a3 



(16) 



Eq. is known as the Dikii-Eilenberger equation<2^ 

Note that we have derived the above equation from 
the BdG equation only with the Andreev approxima- 
tion. On the other hand, it is well known that the Eilen- 
berger equation can be derived from the BdG equation 
or equivalently Gor'kov equation using the quasiclassical 
approximation in addition to the Andreev approxima- 
tion in 3DM In the present case, since we assume the 
system is (quasi) ID where the Andreev approximation 
and the quasiclassical approximation are equivalent, the 
Dikii-Eilenberger equation can be derived without explic- 
itly using the quasiclassical approximation. Therefore, 
the solutions of Eq. (fT6|) are the exact self-consistent so- 
lutions of the BdG equation. 

Next step is to make an ansatz for the form of the 
Gor'kov Green's function. From the gap equation, A(a;) 
must satisfy 
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From Eq. (|10p . the simplest ansatz for R{x] E) to satisfy 
Eq. (|17p is that the diagonal entries are set to be propor- 
tional to I A (a;) p. The gap equation can be derived by 
the functional derivative 



A(a;) cx Tid^e [tTi(l + a3)R{x, E)] 



(18) 



The simplest solution for Eq. is for ofF-diagonal 

entries of the i?(x; E) to be proportional to A{x) (or 
A*(x)). However, the consistency between ([T6l) and (fT8|) 
requires that the derivative term A'{x) (or A*'(x)) should 
be in the off-diagonal entries. The last condition comes 
from Eq. ((T2|) . If we assume the following form for the 
resolvent 
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VYi{a+\A{x)f) bA{x)~icA'{x) 
bA*ix)+icA*'{x) VFt(a-f |A(a;)n 



(19) 



by substituting Eq. ([T9| into the right hand side of 
Eq. (fT6|) . we obtain 



d^R{x; E)a3 = M{E) 



A 



B 

~A 
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icA* 



where 

cdAp)', 

iVFiivp^ +Vpl)E{b{E)A* 
-2ivF^A*{a{E) + lAp). 

Then, we find 

C = UFt^FJ,, 

from the diagonal part, and 

A" + i[b - 2E]A' - 2[a - Eb]A ~ 2A\A\^ = 0, 
from the off diagonal part. Here, we have defined 

a = a^'^a, b = a"^6, E ^ a^'^E, A = a^^A, 
where a is the imbalance parameter, that is 



WFt + "Fi 



(20) 

(21) 
(22) 

(23) 

(24) 

(25) 

(26) 



< a < 1. We introduce the mean Fermi velocity by 
Vf = (vFf + vfi)/2, and set vf — 1- In the balanced 
case (wFt = I'Fj,), it has been shown that this equation 
reproduce the well known solutions, e.g., homogeneous 
condensate (BCS theory), single kink condensate^^, and 
real kink crystal—-. 

To confirm the consistency condition Eq. (I13p . we cal- 
culate the determinant of the resolvent: 



deti? 



./|2 



{2h-P)\A\ 



+z6(A'A* - AA*') 



a2 



(27) 



The NLSE implies that the right hand side of the above 
equation is constant as follows 



(2|A|2 + 2a-62)(|A|2)' 

-(A"A*' + A'A*") 
+ib{A"A* - A*"A) 



0. (28) 



deti? 



It is remarkable that by the scalings in Eq. ([25|) . the 
NLSE for finite spin polarization takes exactly the 
same form as the one for zero spin polarization. This 
means that there exist solutions of Eq. ((24)) which cor- 
respond to each of the solutions of the NLSE for the 
balanced case. Thus, the exact solutions of Eq. ([M)) can 
be easily derived by scaling the solutions of the NLSE for 
the balanced case. The solution corresponding to a com- 
plex kink crystal is the most general one which includes 
other solutions in some limits. We derive the solutions 
of Eq. dMl) in Sec. HVl including that of the LOFF state. 
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III. IMBALANCE EFFECT ON 
SINGLE-PARTICLE STATES 

In this section, we examine the effect of spin imbalance 
on single-particle states. For simplicity, we consider the 
case when the order parameter is real. 



A. Imbalance effect on fermionic zero mode 

First, we consider the fermionic zero mode, i.e. the 
solution of Eq. with E = 0. The zero mode plays 
crucial roles for the LOFF state. The wave function of 
the zero mode is localized around the nodes of the order 
parameter and accommodate the excess spin component. 
For the zero mode solution, we can exactly solve the BdG 
equation. When ii^ = 0, by the scaling transformations 



u = {l + e/2y 



u, 

X = 



ij={l-e/2y^v, (29) 
(1 - e^/4)x = a^x, (30) 



Eq. (U) can be rewritten as 



OX 
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A i# 

L dx -i 




V 



_ 0, (31) 

L ^ ' 95 J L J 

where e = wpf — wpj,- It is clear that the above equation 
has the same form as the one for wpt = wpi- This indi- 
cates that if the BdG equation has a zero mode solution 
for balanced case, there exists a corresponding zero mode 
solution for imbalanced case, and the two solutions are 
related by the scaling transformations (p9)) and (|30l) . 

Furthermore, the zero mode solution can be explicitly 
constructed as follows. If one applies the unitary trans- 
formation 

^ 1+] - — l'^ U " 1 (32) 
Eq. dSI]) yields 

T a] /± = 0. (33) 

Hence, the solution of Eq. (|33)) can be formally written 
as 



f± (x) cx exp 



± / dy a'A{y) 





(34) 



Equation is valid if f± is normalizable. That is 
dxfj. is finite. The solution (p4| is exactly the same 
as the one for the balanced case up to the scaling factor 
for the order parameter. 



B. Perturbation theory for massive modes 

We develop a perturbation theory for massive modes 
{E > 0). We calculate spin imbalance correction for the 



solutions of Eq. (|4]) perturbatively taking e as a small 
parameter. We expand the solution of Eq. (U) by e as 



+ 0(e'),(35) 
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_ w(i)(x) _ 
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i5(0)+gi^(l)+O(,2)_ 



(36) 



In the 0th order, we indeed obtain the equation for 
balanced case 



A(x) 
A(x) ^£ 
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Here, we used A = aA = ^1 -e2/4A w (1 - e^/8)A. 
Making use of the unitary transformation 
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Eq. ([321) becomes 

(-a^ T A' + A2 /f(x) = o. 



(39) 



Once A is obtained by solving Eq. (|24p . the Shrodinger 
type equation ([M]) yields a set of unperturbed eigen- 
states. 

From the first order terms, we obtain 

(-92 ^ A' + A2 - ^(0)) f^^\x) + §a4°)(x) 

+^£;("'9,/i"'(x) - 2£;(")£;(i)/i''^(x) = 0. (40) 

Thus, the first order correction for the energy can be 
calculated from the nonpcrturbative eigenvalues En"^ and 
eigenstates as 



2 



dx 



(0)* 



/l"L(x) 



.(41) 



IV. SPIN IMBALANCE CORRECTION FOR 
VARIOUS CONDENSATES 

A. Homogeneous condensate 

Here we consider the homogeneous condensate 

A(x) = m. (42) 

We can always take m to be real due to the chiral sym- 
metry of the model. The substitution Eq. (|42p into Eq. 
(m yields2^ 



a = 2£;2 - m^, h = 2E. 



(43) 



Next, we calculate the energy spectrum for the quasi- 
particles. In this case, the order parameter is constant 
and thus the zero mode (|34l) is not allowed, i.e. it is 
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not normalizable. Then we calculate the massive modes. 
From Eq. ([37]), we obtain 



''dx 



^(0) 
,(0) 



(x) 



u(0)(a;) 



Then we obtain 

and the eigenspinor is 

,±(0) 
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>{x) 
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where and Vk is independent of x. 
dUl) into Eq. (gf]) yields 



1 



-k. 



(44) 



(45) 



(46) 



Substituting Eq. 



(47) 



Then we obtain the energy dispersion 



£;(fc)± =±VA?2+^+-fc + 0(e2). (48) 

In this case, the spectrum of the BdG equation ^ can 
also be calculated exactly. By substituting Eq. ([25]) and 
Eq. gl]) into Eq. (g]), we obtain 



E±{k) = ±\/k'^ + oP-m?- + -A:. 



(49) 



This dispersion relation is plotted in Fig.[T] When e ^ 1, 
-E ~ ± + TO^, which is consistent with the result 



(j48|) obtained by the perturbation theory. 

When > a^m^, E ~ ±|A:| + efc/2. This indicates 
that the energy dispersion asymptotically becomes that 
of the free fermion (E = wpt^i —vpik). 

Note that the energy gap contracts by a factor com- 
pared to the balanced case, namely the edges of the posi- 
tive and negative energy bands become ±ma^ as plotted 
in Fig. m 



B. Single real kink condensate 

Next solution we consider is the single real kink (anti- 
kink) condensate 



A(x) = ±m tanh(rna;). 
This solution can be obtained by setting 



a = 2E^ - m^, b ^ 2E. 



(50) 



(51) 



Now we analyze the spectrum of the associated BdG 
equation. The energy spectrum for single kink conden- 
sate is obtained in the limit of infinite periodicity of the 
real kink crystal, which will be discussed in the following 




FIG. 1: Fermionic spectrum in the case of homogeneous con- 
densate A = 2 with e = 0.2. The slope of the asymptotes 
changes by e/2 as a consequence of the spin imbalance. The 
dashed line shows the asymtotes for the spectrum when e = 0. 



section. The only exception is that the normalizable zero 
mode exists in this case. For the real kink case (anti-kink 
case), the eigenstate f+{x) {f-{x)) in Eq. ([M]) is normal- 
izable and f{x) {— f+{x) for kink, = ,f-{x) for anti-kink) 
becomes 

f{x) = N[sech{mx)]°'\ (52) 
where N is the normalization constant. 

C. LOFF state 

As shown in Ref. 'Tt!, the NLSE for balanced case has 
the LOFF solution (real kink crystal). Then, we can im- 
mediately conclude that Eq. has the corresponding 
solution 



A(x) 



2m 



isn 



2m 



(53) 



where sn is the Jacobi elliptic function with real elliptic 
parameter < < 1. Then we can conclude that the 
spin imbalance results in the dilatation with a factor 
of the condensate. The substitution Eq. ([53| into Eq. 
p4| yields 



h{E) 
h{E) 



2E^ 
2E. 



2m' 



1 



(54) 
(55) 



The eigenstates of the quasiparticles for the real kink 
crystal order parameter ([55]) are given as foUowsi^, 



pjx + ujs) - 
2L{p-e) 



X exp 



iC (E) 



dx' 



p{x' + uj3)-e 



,(56) 
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where L is the size of the system and p is the Weierstrass 
function which obeys A^(x) — A'(a::) ~ ei + 2p{x + W3) 
with 



2m^ 



ei 



62 



63 



m 

3(1 + 07)2 



r(l + ^), 



3(1 + 



(1 -6^^+^/), 
r(l + 6V^+z^), 



e = ei-£;^ 
The ampUtude of /+^„ has the half periodicity of 

Lu = K{v)/m^ 



(57) 

(58) 

(59) 
(60) 

(61) 



where K{v) is the complete elliptic integral of the first 
kind. We define p as the average of p 



P = - / p{x + ujz)dx. 



The coefiicient C{E) in Eq. ([5E| is defined by 



^(i?) = ±i?^(£;2-i?|)(£;2-£;|), 

= ei-e, (i = 2,3). 
Substituting Eq. into Eq. (liTj) . we obtain 



(62) 



(63) 
(64) 



(65) 



This result implies, in particular, that the two gaps 
shrink as a consequence of non-zero imbalance, as illus- 
trated in Fig. 2. 



E*+eE* 



-E* 
-E2 



-E, 



-E, 



FIG. 2: Schematic illustration of the first order correction to 
the fermionic energy levels with respect to the balanced case 
for the real kink crystal condensate. 



In the limit of — > 1, the periodicity becomes infinite 
and then the LOFF state (the real kink crystal conden- 
sate) reduces to the single kink condensate. 



V. CONCLUSIONS 



We conclude this paper with few remarks. In this pa- 
per, we have investigated the spin imbalance correction 
for the BdG equation. We have expanded the method in 
Ref. Tt*, which is valid for the balanced case, and have 
obtained the non-linear Schrodinger equation for the or- 
der parameter with spin imbalance. We have shown that 
the imbalance correction for the order parameter is only 
included in the reparameterization, however this result is 
nontrivial without using this method. We have obtained 
the fermionic zero-mode exactly for arbitrary spin imbal- 
ance by a scale transformation of the one in the balanced 
case, which implies the stability of the fermionic zero- 
mode against the spin imbalance. We also have analyzed 
the massive fermionic spectrum of the BdG equation with 
small spin imbalance by the perturbation theory. We 
have applied the method for homogeneous condensate, 
the single kink condensate and the LOFF state (the real 
kink crystal condensate). For the homogeneous conden- 
sate, we show the con- sistency between the perturbation 
theory and the exact solution for fermionic spectrum in 
first order. For the real kink crystal condensate, we have 
obtained the imbalance correction for the order parame- 
ter and the fermionic spectrum. This result completely 
generalizes those of Ref. 10, 

Finally, we make few remarks, (i) we have analyzed the 
fermionic problem by the perturbation at the first order 
in the spin imbalance parameter e, however the higher 
expansion is straightforward, (ii) As in the case of ho- 
mogeneous condensate, we may obtain the exact solution, 
(iii) We have dealt with the real condensate in this paper, 
however the method used here can be generalized to the 
complex case, such as the twisted kink^^ and the twisted 
kink crystalii. 
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